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Abstract. We present a consistent model for the UV 
and supersoft X-ray emission from the symbiotic nova 
SMC3 (= RX J0048.4-7332). Following the present pic- 
ture of symbiotic stars, the model consists of radiation 
from a hot star and an emission nebula excited by that 
star. The observations were compared to theoretical mod- 
els in which the hot star's emission is calculated with the 
help of hydrostatic and Wolf-Rayet-type non-LTE model 
atmospheres. Our analysis clearly shows evidence for mass 
loss rates of several 10 -6 M /yr. The minimum effective 
temperature compatible with both the observed UV and 
X-ray flux is about 260 000 K, which is higher than in any 
other star analyzed with sophisticated NLTE model at- 
mospheres. Since the hydrostatic surface is hidden by the 
stellar wind no upper limit for the temperature can be 
determined. However, we were able to determine the to- 
tal luminosity of a symbiotic nova with reasonable accu- 
racy (LsMC3 = 1O 4 O5±O O5 L ). This value is well below 
the Eddington limit (« 50 000 L ). In order to reproduce 
the observed energy distribution a carbon-to-hclium ratio 
> 2 ■ 10~ 4 — leading to an absorption edge at 0.39 keV — 
is necessary. 

Key words: Binaries: symbiotic - Stars: fundamental 
parameters - Stars: Individual: SMC3 - Stars: Individual: 
RX J0048. 4-7332 - Stars: mass-loss - X-rays: stars 



1. Introduction 

ROSAT observations established the class of 'supersoft X- 
ray sources' ('SSS'; e.g. Hasinger 1994) with almost no flux 
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* Based on data collected with the ROSAT and the IUE 
satellites 



at hv > 0.5 keV; normally, close-by low-luminosity objects 
like single white dwarfs or cool coronal stars are excluded 
from this definition. The nature of many SSS is still a mat- 
ter of controversial discussion. RX J0048. 4-7332 is one of 
ten SSS confirmed as members of the Magellanic Clouds 
(Kahabka & Pietsch 1993; Pietsch & Kahabka 1993). In 
this paper we investigate the hypothesis that the X-ray 
flux of this object is due to photosphcric emission of a 
very hot white dwarf. Given a sufficiently high temper- 
ature the photospheric emission of white dwarfs is mea- 
surable in the ROSAT window (e.g. Jordan et al. 1994b, 
Wolff et al. 1996). A hot white dwarf is hardly the cor- 
rect explanation for all supersoft X-ray sources. However, 
at least for some sources, the presence of a very hot star 
is also in agreement with other evidence: For instance, in 
symbiotic systems there are high excitation nebular emis- 
sion lines combined with a hot continuum in the UV. The 
supersoft appearance is also expected from the theoretical 
side; in particular, for symbiotic novae Sion & Starrficld 
(1994) predict phases in which a very hot white dwarf can 
be detected. 

Indeed, the supersoft source RX J0048. 4-7332 coin- 
cides with the symbiotic nova SMC3 that has been in 
outburst since 1981 (Morgan 1992). Thus, for this sys- 
tem it appears natural to propose that a hot white dwarf 
is the source of the X-ray flux. Symbiotic stars are inter- 
acting binaries consisting of a red giant and a very hot 
white dwarf (typically T e g ~ 100 000 K), whose radiation 
ionizes the emission nebula. In the optical spectral range 
we observe a composite spectrum from the cool star and 
the nebula. The UV spectrum is dominated by the nebular 
emission with only a small contribution from the hot star 
that becomes dominant at the shortest wavelengths. The 
spectral range best suited for a direct observation of the 
hot component of a symbiotic system is the soft X-ray and 
EUV part of the electromagnetic spectrum, where these 
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stars emit the bulk of their energy (e.g. RR Tel, see Jordan 
et al. 1994a, hereafter JMW). 

Vogel & Morgan (1994; Paper I) obtained an IUE spec- 
trum of SMC3 which looks typical for a symbiotic star 
of moderate excitation. On the other hand, the optical 
spectrum (Morgan 1992; Miirset et al. 1995 = Paper III) 
reveals nebular ionization that is unusually high for sym- 
biotics, up to Fe +9 . Even more peculiar, however, are the 
X-ray observations presented by Kahabka et al. (1994): 
Despite of the large distance, the ROSAT count rate of 
SMC3 is about the same as of RR Tel, hitherto the X-ray 
brightest of all known galactic symbiotic novae (Jordan et 
al. 1994a; hereafter JMW). Thus, the ratio of the optical 
to the X-ray flux of SMC3 is a challenge for models of 
symbiotic stars. On the other hand, the known distance 
of this extragalactic system provides a unique opportunity 
to investigate the energetics of interactions and outbursts 
in symbiotic systems. 

Kahabka et al. (1994) have presented black body fits 
to the ROSAT data, which resulted in an extremely unre- 
alistic luminosity of up to 10 13 L Q . However, we will show 
that a reasonable value is obtained if the data are com- 
pared to appropriate model atmospheres, and if we con- 
sistently take into account the constraints from both the 
X-ray and UV regions of the electromagnetic spectrum. 

In Sect, p we list the data used in our analysis, in Sec- 
tions U and | we describe the methods, the results are 
discussed in Sect. M, and Sect. || presents some concluding 
remarks. 



2. Data 

2.1. ROSAT 

Pointed ROSAT PSPC observations were retrieved from 
the ROSAT archive and reduced using the EXSAS soft- 
ware (Zimmermann et al. 1994). Kahabka et al. (1994) 
describe the data, obtained during a 20 ks exposure in 
April 1992, in detail. We binned the counts into the 
various energy channels until a signal-to-noise ratio of 
five was achieved. After correcting for dead time and vi- 
gnetting the resulting total count rate amounts to 0.204 ± 
0.004 counts/sec. The observed PSPC pulse height distri- 
bution (PHD) is displayed in Figs. || - [| 

A follow-up observation in December 1992 is not in- 
cluded in the archive by the time of writing, and could 
therefore not be used for this work. 

2.2. IUE 

The UV spectrum of SMC3 is described in Paper I. It 
was observed with the short wavelength (1200 A A 5s 
1900 A) camera of the International Ultraviolet Explorer 
satellite (IUE) at low resolution using the large (20" x 10") 
aperture in April 1993. Given the stability of SMC3 in 
its present outburst stage (Morgan 1992) we assume that 



there was no variability so that the ROSAT and the UV 
observation can be combined. 



2.3. Optical spectra 

Optical low resolution spectra of SMC3 have been pub- 
lished by Morgan (1992) and in Paper III. The spectrum 
is typical for a high excitation s-type symbiotic. It displays 
Hi Balmer, He I, and Hen emission lines, as well as the 
Raman scattered Ovi A A 6825, 7082 features (cf. Schmid 
1989). In addition to the typical features, [Fex] A6374 
is also present which is very rarely observed in symbiotic 
spectra. Davidsen et al. (1977) report [Fex] emission from 
GX1+4, which is a symbiotic system that contains a neu- 
tron star, and Webster & Allen (1975) observed [Fex] for 
He2-171. The spectrum of AS 295B contains [Fex], [Fexi], 
and [Fexiv] (Herbig & Hoffleit 1975). 

From the H/3 luminosity and the line ratios given by 
Morgan (1992) we derive the optical line fluxes of H/3 and 
He II A4686 listed in Table 0. 



Table 1. Continuum and integrated emission line fluxes from 
SMC 3 used in the analysis. The optical line fluxes are from 
Morgan (1992). 



/(continuum A = 1300 A) 2.9 ■ 10 

/(HenA1640) 2.3-10 
/(Hen A4686) 
/(H/3 A4861) 

/([Fex] A6374) 2.0 • 10~ 14 erg cm' 2 s" 1 



" 15 erg cm 2 s 1 A 1 

. i ■a _o _ 1 

erg cm s 
3.3 • 10 -14 erg cm" 2 s" 1 
4.5 • 10~ 14 erg cm" 2 s" 1 



2.4- Distance and interstellar extinction 

Our analysis requires reddening corrected continuum and 
emission line fluxes. The extinction can be derived by com- 
paring the flux ratio of the He n recombination lines A4686 
and A1640 to the theoretical value for case-B recombina- 
tion (Hummer & Storey 1987). Although this ratio im- 
plies negligible extinction, we corrected for the foreground 
extinction of the SMC {E B -v = 0™05; Bessell 1991) by 
using a galactic extinction curve, and subsequently apply- 
ing the SMC extinction curve by Prevot et al. (1984) for 
E'b-v = 0^03, which is the average internal extinction of 
the SMC (Bessell 1991). For the distance to the SMC we 
adopt d = 61 kpc determined by Barnes et al. (1993). 

3. Atmospheres 

The spectroscopic analysis is based on state-of-the-art 
non-LTE atmospheres as required by the high tempera- 
tures involved. They were calculated with two different 
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codes, one assuming plane-parallel geometry and hydro- 
static layers, the other is appropriate for extended atmo- 
spheres with a steady state stellar wind. Table ^ gives an 
overview on the hot model atmospheres that have been 
calculated. 

3.1. Hydrostatic models 

As in the case of RR Tel (JMW) we first computed the 
theoretical model spectra for the hot stellar component 
with line blanketed non-LTE model atmospheres assum- 
ing plane-parallel geometry and hydrostatic and radiative 
equilibrium. Besides hydrogen and helium the models con- 
tain C, N, and O. The number of non-LTE levels assigned 
to each species are 5, 12, 38, 48, 48, respectively, and a 
total of 333 lines were explicitly included. The character- 
istics of these models were described in detail by Werner 
& Heber (1991). 

3.2. Models with a stellar wind 

Plane-parallel models are justified if the thickness of the 
stellar atmosphere is small compared to the stellar radius. 
This is normally the case for hot white dwarfs. However, 
these assumptions may not be valid in the case of a sym- 
biotic nova, since for several of the symbiotic novae, there 
is evidence for mass loss (e.g. Nussbaumer et al. 1995, 
Nussbaumer & Vogel 1995). Therefore, we also calculated 
models with mass loss as an additional free parameter us- 
ing the co-moving frame code for hot stars with a strong 
stellar wind developed in Kiel (Hamann & Schmutz 1987; 
Wessolowski et al. 1988). The non-LTE calculations in- 
clude 6 levels for the hydrogen model atom; 11 levels for 
helium with 5 levels for He I and 5 levels for Hen; 12 lev- 
els for carbon with 5 levels for C iv and 6 levels for C v; 6 
levels for nitrogen (5 levels for N v) and 9 levels for oxygen 
(5 levels for O vi). In total 44 levels and 54 line transitions 
are treated explicitly. 

The free parameter of a model with an extended at- 
mosphere due to mass loss are T e g, i?*, M /foo, and (3. 
The radius, i£*, is the inner boundary of the model. This 
inner boundary is either where the expansion velocity of 
the envelope is subsonic or, if the mass loss is so strong 
that this region is at very large optical depths, then the 
inner boundary is set where the Rosseland optical depth 
tr « 100. The temperature T c g yields in combination with 
i?* the luminosity. Voo and (3 specify the velocity law of the 
wind. We adopt a velocity law v(r) = Uqo(1 — R*/ r ) 13 , as- 
suming (3 = 1 and = 1000 km/sec (cf. section 5.2). This 
type of law resembles that found empirically for O stars 
(Groenewegen & Lamers 1989) and agrees with the law 
calculated hydrodynamically for radiation driven winds of 
O stars (Pauldrach et al. 1986; Schaerer & Schmutz 1995). 



4. Fit procedure 

There are three observational constraints that have to be 
reproduced by a successful model: the X-ray flux, the UV 
continuum in the IUE wavelength range, and the nebular 
emission line Hen A1640; in addition, there is a high ion- 
ization line, [Fex] A6374. Since a forbidden transition is 
very sensitive to the unknown density of the nebula, the 
observation of the [Fe x] line implies only a lower limit for 
the hardness of the exciting radiation, but the observed 
line strength does not help to constrain a stellar parame- 
ter. 

In order to reproduce the three observations we have 
adapted the following procedure. First, we have computed 
a large grid of model atmospheres with different param- 
eters in T c ff, logy, chemical composition, and mass loss 
rates (or M /woo, respectively). 

The emergent flux does not depend on the radius of 
the models. This is obviously true for hydrostatic models 
where the radius is not among the model parameters. But, 
it is also correct for spherically extended models with mass 
loss, provided that the product M /voo/R^ 2 is kept con- 
stant (Schmutz ct al. 1989). Thus, if we assume that the 
continuum at A = 1300 A is dominated by the flux of the 
hot star, we can compute the stellar radius for each model, 
knowing the flux predicted by the model atmosphere and 
the observed luminosity at 1300 A; at this wavelength the 
hot stellar emission dominates the continuum in all sym- 
biotic objects whose UV spectra have been investigated 
by Miirset et al. (1991). 

With the absolute dimensions of the star determined 
by the observed 1300 A flux for each model of the grid, we 
fed the emergent radiation to the photo-ionization code 
NEBEL (Nussbaumer & Schild 1981; Deuel 1986; Vogel 
1990). For the calculation of the emission lines we assume 
that the nebula is radiation bounded in the Hen region. 
We have no possibility to verify this assumption. However, 
Miirset et al. (1991) found the Hen region to be radiation 
bounded for almost all symbiotic systems and only a few 
exceptions being nearly radiation bounded. The presence 
of a radiation bounded nebula is in no contradiction to the 
observed X-rays since our calculations show that an ion- 
ization bounded nebula is not necessarily optically thick 
below 100 A; material from the neutral part of the nebula 
does, of course, contribute to the column density towards 
the X-ray source. Indeed, we find hydrogen column den- 
sities (up to 1.1 ■ 10 21 cm -2 ) that are much larger than 
the value we would expect from the reddening of SMC3 
(the formula by Groenewegen & Lamers 1989 yields about 
iV H = 1.9- 10 20 cm' 2 for £ B -v = m 05). We interpret the 
large hydrogen column densities as evidence for the sym- 
biotic nebula being radiation bounded. 

The strength of the He n A1640 line is a measure for 
the number of stellar photons emitted below A = 228 A, 
but weighted strongly for photons just shortwards of the 
threshold wavelength. Thus, basically, we use a sophisti- 
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cated Zanstra method. As it is the case for the original 
Zanstra method, our analysis depends only weekly on the 
adopted parameters of the nebula. A comparison of the 
calculated He II A1640 line flux with the observed one se- 
lects those models with the correct number of He + ionizing 
photons. 

For a black body energy distribution the Hen A1640 
analysis yields a unique temperature. If we use emergent 
radiation from hydrostatic atmospheres, the Hen A1640 
analysis selects different effective temperatures depending 
on the chemical composition of the atmosphere and - to 
a much lower degree - on the gravity. However, although 
there is a range of temperatures of plane-parallel models 
that fits the observed Hen A1640 flux, the Hen A1640 
analysis still determines the stellar temperature reason- 
ably well. 

The results using atmospheres with a stellar wind differ 
fundamentally from those based on plane-parallel models. 
If we feed the model flux resulting from our calculations 
with extended model atmospheres with mass loss to the 
He II A1640 analysis, we find that the number of He + ioniz- 
ing photons is a sensitive function mainly of the mass loss, 
and not of the stellar temperature (see Sect. 5.3). The flux 
below 228 A is suppressed by absorption in the wind de- 
pending on the mass loss rate. If the mass loss rate is suf- 
ficiently large helium recombines to He + and in this case, 
the wind is optically thick to photons more energetic than 
the He+ edge (Schmutz et al. 1992). In Fig. | we show 
the EUV and soft X-ray energy distribution predicted by 
three models with about the same effective temperature 
of T c ff = 300, 000 K but different mass loss rates (models x 
- z). Thus, for wind models the Hen A1640 analysis does 
not provide information on the effective temperature ex- 
cept that there is a lower limit in order to produce He + 
ionizing radiation. However, no upper limit to the effec- 
tive temperature of the star can be specified: We can al- 
ways increase the mass loss rate to reduce the number of 
He + ionizing photons to the value implied by the observed 
He II A1640 flux, no matter what stellar temperature we 
choose. 

The third step of our analysis is to compare the X- 
ray flux of the model atmosphere with the flux observed 
by ROSAT. For this comparison the theoretical spectrum 
has to be attenuated to account for the absorbing inter- 
stellar and circumstellar material. We calculate the ab- 
sorption according to Morrison & McCammon (1983) and 
Cruddace et al. (1974). The attenuated spectrum is then 
converted to a pulse height distribution (PHD) by convo- 
lution with the PSPC detector response matrix and the 
effective areas provided by the EXSAS package. The com- 
parison of the theoretical PHD with the low energy part 
of the ROSAT observation, 0.1 to 0.2 keV, allows to de- 
termine the hydrogen column density. The comparison of 
the high energy part tests the emergent flux of our mod- 
els at wavelengths where the flux is sensitive to the stellar 




100 200 
wavelength / A 



Fig. 1. Relative energy distributions of three models 
with approximately the same effective temperature of 
T e ff= 300 kK but with different mass loss rates (model 
x: M £ = 10" 5 1 M /yr, y: My_= 10~ 5 ' 2 M /yr, and z: 

M z = 10~ 5 ' 3 M /yr; see Table|| for the other parameters). 
Note that a final velocity of Voo = 1000 km/sec is assumed 
(only M /vca enters as a parameter into the model calcula- 
tion) 



temperature and to the presence of metals (e.g. the C +4 
ionization edge at 0.39 keV). 

5. Results 

In Table [2] we list those parameter ranges for which we 
have calculated model atmospheres and for which the 
prediced flux distribution is able to reproduce the observed 
continuum flux at 1300 A as well as the Hen A1640 line 
flux. For these models the predicted X-ray luminosity was 
compared to the ROSAT observations. In Table ||we sum- 
marize the results of these fits, and list the resulting T ff 
and £*, and for the models with stellar wind M /voo- Also 
shown is the hydrogen column density required to fit the 
low energy X-ray part. In principle we could quantify the 
quality of the fits by listing the reduced x 2 of the solution. 
However, one can clearly see that the deviations between 
the predictions and the observations are not solely due to 
statistical errors: Even in the case of the "relative good 
fits" systematic deviations remain which may have its ori- 
gin in small incompleteness of the modelling. Finally, Ta- 
ble H comments on the quality of the resulting fit to the 
"high" -energy tail X-ray flux. 

5.1. Black body emission 

If we assume that the radiation of the hot component has a 
black body spectral distribution the Hen A1640 analysis 
yields a unique solution, namely Tbb = 100, 000 K and 
R = 0.3 R©. However, the emergent X-ray flux of this 
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Table 2. Summary of model atmospheres parameters. In each series we computed atmospheres of several temperatures until a 
fit to the UV continuum and the He II emission could be achieved. The chemical abundances are given as number fractions. 



Series 


Model type 


log M /Voo 

[M /(yr km sec" 1 )] 


log 


9 




H 


He 


C 


N 


O 


Abundances 
















[io- 5 ] 


[IO" 5 ] 


[IO" 5 ] 


similar to 


A 


hydrostatic 


- 


5 - 


7 


- 


0.90 


0.10 


48.00 


10.00 


84 


Sun 


B 


hydrostatic 


- 


6 




- 


0.90 


0.10 


4.80 


1.00 


8.4 


Magellanic Clouds 


C 


hydrostatic 


- 


7.5 




- 


1.00 


0.00 


0.00 


0.00 


0.00 


DA 


D 


hydrostatic 


- 


6 - 


7.5 




0.01 


0.99 


100.00 


0.00 


0.00 


DO 


E 


hydrostatic 


- 


6 - 


7 


- 


0.01 


0.99 


0.10 


0.00 


0.00 


DO 


F 


hydrostatic 


- 


5.5 


- 6 


- 


0.67 


0.33 


10.00 


0.00 


0.00 


DO 


G 


hydrostatic 


- 


6 






0.50 


0.50 


100.00 


0.00 


0.00 


DO 


H 


hydrostatic 


- 


6 - 


7 


- 


0.50 


0.50 


10.00 


0.00 


0.00 


DO 


P 


WR-type 


-8.53 


~ 




1 


0.50 


0.50 


40.00 


0.76 


7.50 


~ 


q 


WK-type 


—8.50 






1 


0.50 


0.50 


40.00 


0.76 


7.50 




r 


WR-type 


-8.12 






1 


0.50 


0.50 


100.00 


0.76 


100.00 




s 


WR-type 


-7.97 






1 


0.50 


0.50 


200.00 


0.76 


100.00 




t 


WR-type 


-8.30 






1 


0.50 


0.50 


20.00 


0.76 


7.50 




u 


WR-type 


-8.30 






1 


0.50 


0.50 


40.00 


0.76 


7.50 




X 


WR-type 


-8.10 






1 


0.50 


0.50 


40.00 


0.76 


7.50 




y 


WR-type 


-8.20 






1 


0.50 


0.50 


40.00 


0.76 


7.50 




z 


WR-type 


-8.30 






1 


0.50 


0.50 


40.00 


0.76 


7.50 





Table 3. Solutions to the UV continuum and Hell A1640 emission. The 5' column tells us whether a fit to the ROSAT 
observation could be achieved. The comment 'too soft' means that too much counts are predicted at low energies compared to 
the observation. Column 6 lists the implied interstellar hydrogen column density. 



Model 


Model type 


Teft 




Comment on implied ROSAT solution 


Ah Fig. 






[kK] 


[L©] 




[10 20 cm- 2 ] 


Black Body 


Black Body 


100 


8100 


count rate is by far too low 


< 0.1 



A 
B 
C 
D 
E 
F 
G 
H 



hydrostatic 100 11000 count rate is too low 

hydrostatic 120 17000 synthetic PHD is too soft 

hydrostatic 65 5 800 synthetic PHD too low at 0.3-0.4 keV 

hydrostatic 130 23 000 synthetic PHD too low at E J> 0.25 keV 

hydrostatic 120 17000 synthetic PHD too low at E > 0.20 keV 

hydrostatic 120 17000 synthetic PHD too low at E Z 0.20 keV 

hydrostatic 120 16 000 synthetic PHD too low at E <; 0.20 keV 

hydrostatic 120 17000 synthetic PHD too low at E £ 0.20 keV 



WR-type 220 9 200 PHD slightly too low at 0.20-0.45 keV 

WR-type 265 10 200 PHD slightly too low at 0.30-0.45 keV 

WR-type 309 11200 relative good fit 

WR-type 354 12 700 relative good fit 

WR-type 300 16 800 PHD too high at E > 0.45 keV 

WR-type 300 16 800 PHD slightly too low at 0.30-0.42 keV 

WR-type 315 11500 PHD too high at E > 0.37 keV 

WR-type 325 16 000 PHD too high at E > 0.35 keV 

WR-type 335 22 900 PHD too high at E > 0.33 



0.0 
0.1 
0.4 
4.0 
1.7 
2.0 
2.1 
2.6 



P 

q 

r 

s 
t 
u 

X 

y 

z 



7.7 
8.7 
8.8 
9.3 
8.8 
8.8 
11.0 
11.0 
11.0 



solution is much too low to reproduce the PSPC count 
rate, even for a vanishing interstellar hydrogen column 
density. Moreover, no [Fe x] emission is expected because 
the radiation field is not hard enough to produce Fe +9 . 

Good fits to the PHD can be achieved if cither the stel- 
lar radius is increased to about 1.5 • 10 7 Kq or the temper- 
ature is increased to T BB £ 220,000K (with R~2 R Q ), 



in agreement with the result seen in Fig. 1 of Kahabka et 
al. (1994). On obvious reasons, these parameters clearly 
lead to an U V flux several orders of magnitude too high 
compared to the observation. Also, the luminosity would 
exceed the Eddington limit by many orders of magnitude. 
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0-30 0.40 
channel energy / keV 



Fig. 2. The solid line is the best solution to the UV and 
X-ray observations with DO model atmospheres (model D 
in Tab.|): He/H = 100, C/H = 0.1, T cff = 130, 000 K, and 
Nn — 4.0 • 10 20 cm -2 . Also shown is one example for another 
chemical composition (dashed line, model E): He/H = 100, 
C/H = 10 -4 , T cff = 120, 000 K, and iV H = 1.7- 10 20 cm" 2 . The 
other compositions result in nearly the same synthetic PHDs 



5.2. Hydrostatic models 

Models with solar abundances (series A, see Tab. 3) lead 
to results similar to the black body case: the UV measure- 
ments would result in T e g w 100, 000 K whereas T e g ^ 
200, 000 K is needed to reproduce the X-ray observations. 
If we reduce the CNO abundances by a factor 10 (series 
B), the UV measurements are in agreement with the pre- 
dictions at T c g k 120,000 K. However, this is still too cold 
to account for the X-ray flux. 

Pure hydrogen (DA, series C) atmospheres produce by 
far the "hardest" possible spectrum, i.e. the highest X-ray 
flux for a given temperature. Again, these models lead to 
the situation that a relatively good fit to the PSPC PHD 
can be obtained at T cS » 100,000K, but the Hen A1640 
line flux is only reproduced at T e g <J 65, 000 K. 

Next, we assumed that helium is the main con- 
stituent of the hot component's atmosphere. We exam- 
ined DO/DOA model atmospheres with different helium- 
to-hydrogen ratios. The best solution to both, the UV 
and X-ray measurements is obtained with Hc/H = 100, 
C/H = 0.1, and T cff = 130, 000 K (Fig.|, model D). If 
the hydrogen column density is chosen so that the lowest 
channels of the pulse height distribution are reproduced 
then the predicted flux is somewhat too small at higher 
energies. This deviation is larger for all other combinations 
of helium and carbon abundances (see Fig.[|). In the at- 
mosphere of a hot DO white dwarf with T g > 100, 000 K 
nitrogen, oxygen, and iron group elements may also be 
present. However, if we assume metal abundances as high 



as in PG 1034+001 (T off « 100, 000 K, Werner et al. 
1995a) the X-ray flux is completely absorbed. 
Therefore an ionized spectrum is required 

— with a strong Hen A228 absorption edge in order to 
prevent moderately ionized lines like Hen A 1640 from 
becoming too strong. 

— and a considerable portion of the energy released be- 
yond the "gap" due to Hen and possibly metals, to 
absorb high-energy photons. 

5.3. Models with a stellar wind 

If we assume that the hot object in SMC3 has a stellar 
wind we have to determine three model parameters, one 
more - the mass loss rate - in addition to effective tem- 
perature and stellar radius. As pointed out in Sect. [|, the 
stellar wind allows to suppress the He + ionizing photons 
even for very high stellar temperatures. With three model 
parameters — T e g, i?*, and M /voo — and three obser- 
vational constraints — the continuum flux at 1300 A, the 
flux of the He n A1640 line, and the X-ray emission — it 
should be possible, in principle, to determine a solution. 

However, in practice, we find that no unique solution 
exists. The problem is that if we fit both UV constraints, 
the 1300 A continuum and the Hen A1640 line flux by 
adjusting two of the three parameters, say and M 
/voo as a function of T c g, we find that the predicted X- 
ray flux depends on the third parameter, T c g, only in a 
limited range of temperatures: In particular, we find that 
for lower temperatures (r* 100 — 200 kK) the X-ray flux is 
quite sensitive to T e g; at higher temperatures, however, a 
saturation effect occurs and we obtain an upper limit to 
the predicted X-ray luminosity, almost regardless of the 
specified effective temperature. 

It turns out that this upper limit reproduces almost, 
but not perfectly, the observed X-ray spectrum. The rea- 
son for this result is a typical property of a stellar wind. 
The wind that suppresses the photons at the the He + ion- 
ization edge is still relatively opaque in the observed X-ray 
region below 100 A. This is simply because the bound-free 
continuum cross section is proportional to A 3 . Therefore, 
as long as we have a stellar temperature where the maxi- 
mum of the radiation is in the observed X-ray or beyond, 
we find that all models have similar X-ray energy distri- 
bution determined by the helium absorption cross section. 
This is similar to adjusting the interstellar column den- 
sity Nn. All these arguments are only valid for the energy 
spectrum around 100 A (where the bulk of the observed 
X-ray flux originates). The high energy tail, where the 
influence on the helium bound-free absorption becomes 
smaller, would in principle be very sensitive to the temper- 
ature. However, as discussed below, the presence of carbon 
absorbs all flux below 31.6 A. 

Our difficulty to determine a unique set of the three 
parameters T e s, i?*, and dM/dt does not prevent us to de- 
termine the stellar luminosity. We have stated above that 
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the form of the energy distribution is well determined and 
therefore, since we know the distance to our object, we 
know the luminosity with only a weak dependence on the 
other parameters. The statement that the effective tem- 
perature is not well known, is physically speaking a limit 
in determining the hydrostatic radius of the hot object, 
because this radius is hidden by the stellar wind. From 
models r and s (see Table ||, Fig. ||) that yield the best fit 
to the observations, we find ^SMC3 = 10 4 ' 05±0 ' 05 L Q . 




0.30 0.40 
channel energy / keV 



Fig. 3. Comparison of the observed PHD of SMC3 with three 
Wolf-Rayet type models with T e g — 220 kK (model p, see Ta- 
ble! for the other parameters), 309 kK (r), and 354 kK (s). The 
predicted flux from model p is not hard enough to reproduce 
the observation, while the two hotter models (r) and (s) lead 
to relatively good fit to the data 




0.40 

energy / keV 



Fig. 4. The "coolest" Wolf-Rayet type model with T cff =260 K 
(q) that is able to reproduce the observed ROSAT PHD rea- 
sonably well 
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Fig. 5. Predicted ROSAT PHD for the models with differ- 
ent mass-loss rates 
shown in Fig.|: M E = 10" 5 1 M Q /yr, M y = 10" 5 ' 2 M Q /yr, and 

M z = 10~ 5 ' 3 M Q /yr; Voo = 1000 km/sec is assumed 




0.30 0.40 
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Fig. 6. These wind models with T eS » 270 kK (model t and u) 
show, that too much hard X-ray photons are predicted if the 
carbon abundance is too low (C/He=2 • 10~ 4 in the case of t 
compared to and 4 • 10~ 4 for u) 



Although we cannot decide which temperature is to be 
preferred, T c g = 309 kK, 354 kK, or even hotter, we can 
set a lower limit to the effective temperature, i.e. an up- 
per limit to the stellar radius. In Fig. || we show fits to the 
observed X-ray data based on models with effective tem- 
peratures of T cff = 220 kK (model _p), 309 kK (r), and 
354 kK (s) . ft is evident that the radiation of the model p 
is not quite hard enough to reproduce the observed pulse 
height distribution. On the other hand, both hotter mod- 
els come close to fitting the observed X-rays. Fig. || shows 
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that a model with T c s — 265 kK (model q) is the coolest 
to fit the energy distribution of the hot object in SMC3 
reasonably well. 

In order to derive M a value for d m is required. Unfor- 
tunately, Woo is not easily measurable. We do not dispose 
of the high S/N, high resolution UV spectra which al- 
low to detect broad line components and determine their 
width. Even for galactic objects the velocity of the winds 
from hot components of symbiotic systems could be mea- 
sured only for a few objects. Nussbaumer et al. (1995) 
and Nussbaumer & Vogel (1995) derive Voc ~ 1000 km/s 
which we will adopt in the following for SMC3. With this 
value, model q implies a rather high mass loss rate of 
M> 10~ 5 - 3 M yr- 1 . 

Helium and carbon are the only elements that have 
a significant influence on the resulting spectrum. Helium 
is needed in order to produce the absorption at the He + 
edge. The mass loss rate quoted above is basically the mass 
loss of helium. We have assumed a hydrogen abundance 
of 20 % by mass, but the hydrogen abundance could be 
changed arbitrarily without observable consequences. The 
carbon abundance is found to be very important. With 
the hot temperatures of our best fitting models we predict 
significant flux shortwards of the C +4 ionization edge at 
31.6 A (0.39 keV). Such energy distributions would pro- 
duce much larger count rates than observed at energies 
larger than 0.4 keV. Therefore a relatively high carbon 
abundance is required in order to fit the observed energy 
distribution. Basically, we are invoking the same wind ab- 
sorption effect at the C +4 edge as we invoked above to re- 
duce the numbers of He + ionizing photons. Fig.|^ demon- 
strates the importance of the carbon edge by comparing 
the predicted ROSAT PHD for C/He=2 • 10~ 4 (too much 
"hard" X-ray flux) and 4 • 10~ 4 (slightly too low flux at 
E > 0.3 keV) for about the same effective temperature 
(« 270 kK). A slightly larger relative carbon abundance 
is required to suppress the high energy photons for the 
hotter models. Therefore this carbon abundance is a min- 
imum value and we cannot exclude higher values. 

6. Conclusions 

SMC3 is extraordinary, even among symbiotic stars. It is 
remarkably X-ray bright and shows unusually high ioniza- 
tion in the nebula. We succeeded to model these proper- 
ties and the UV flux with a very hot WR-type model at- 
mosphere. Our result clearly confirms that at least some 
of the supersoft X-ray sources contain very hot photo- 
spheres of degenerate stars. It shows in particular, that 
such hot stars can be encountered in symbiotic binary sys- 
tems. Reasonable fits are achieved at temperatures above 
260 000 K and mass-loss rates of M £1 ■ KT 6 M /yr, (as- 
suming Woo = 1000 km/sec). Due to the properties of the 
wind models (the hot star itself is hidden by the stellar 
wind) no upper limit for the effective temperature could 



be specified. However, this uncertainty does not effect the 
total luminosity: ^SMC3 = 10 4 05±0 05 L Q . 

Even the lower limit for the effective temperature is 
extraordinarily high (see the temperature distribution for 
hot components of symbiotic systems in Fig. 7 of Pa- 
per III). To our knowledge SMC3 has the hottest stellar 
atmosphere ever derived in a sophisticated NLTE anal- 
ysis. Until now, the hottest known (pre-) white dwarfs 
are the DA central star of the planetary nebula WDHS1 
with T off =160. . .200kK (Liebert et al. 1994; Napiwotzki 
1995) and the PG1159 stars H1504+65 (Werner 1991) and 
RXJ2117+3412 (Werner et al. 1995b) with T cff =170kK. 
Note, however, that the temperature of SMC3 refers to 
the hydrostatic radius i?* and not to tr = 2/3. 

On the basis of LTE model atmospheres Heise et al. 
(1994) have shown that the ROSAT PSPC observations 
of the super-soft X-ray source 1E0056. 8-7154, associated 
with the planetary nebula N67 (Wang 1991; Cowley et al. 
1995) can be explained assuming a hot pre-white dwarf 
with an effective temperature as high as 450 kK. However, 
at such high temperatures NLTE effects are expected to 
be strong so that this result has to be checked using NLTE 
model atmospheres. 

Strong mass loss is indispensible to explain the spec- 
trum of SMC3. Direct or indirect signatures of strong mass 
loss are also encountered in other symbiotic novae (see 
Sion et al. 1993, JMW, Nussbaumer et al. 1995, Nuss- 
baumer & Vogel 1995, and the discussion in Miirset et al. 
1995a). If SMC3 has been shedding its wind constantly 
since the beginning of the outburst, it must have lost at 
least a mass of AM ~ 10~ 4 M Q . This is much larger 
than the mass burnt during the same time, which means 
that the duration of the outburst of SMC3 will be limited 
rather by the mass loss than by the nuclear processes. AM 
is also a lower limit to the mass of the layer accreted by 
the white dwarf prior to outburst. 

Kenyon et al. (1993), Miirset & Nussbaumer (1994), 
and JMW found for the known galactic symbiotic novae 
luminosities clearly below the Eddington limit. This result 
is now confirmed with a result that, due to the well known 
distance, is probably more reliable than the galactic ones. 
Here we encounter a basic difference between symbiotic 
and classical novae which seem to exceed the Eddington 
limit (Starrfield et al. 1993; Shore, personal communica- 
tion). If a core- mass - luminosity relation holds for the 
symbiotic novae we can derive a mass M ss 0.8 M Q for 
the hot star (with the relation from Joss et al. 1987). 

As a final comment we would like to stress that a pho- 
tospheric X-ray source like SMC3 can only be reasonably 
investigated with i) multi-frequency data and ii) sophisti- 
cated non-LTE atmosphere models. Kahabka et al. (1994) 
impressively demonstrated the failure of black body fits to 
the ROSAT PHD by deriving a luminosity of L ~ 10 13 L 
for SMC 3. The reasons for the failure are obvious: Fit- 
ting the ROSAT data alone is a dangerous extrapolation, 
because they represent only an extreme tail of the stel- 
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lar spectrum. Even worse, this extreme part is exactly the 
spectral region where hot atmospheres differ very strongly 
from Planck curves. 
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